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properties that have been observed in the vicinity of the transition. 
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1. INTRODUCTION 

In two-dimensional electron systems, the electrons are confined to move in a plane in 
the presence of a random potential. According to the scaling theory of localization 
(Abrahams et al 1979), these systems lie on the boundary between high and low 
dimensions insofar as the metal-insulator transition is concerned. The carriers are 
always strongly localized in one dimension, while in three dimensions, the electronic 
states can be either localized or extended. In the case of two dimensions the 
electrons may conduct well at room temperature, but a weak logarithmic increase 
of the resistance is expected as the temperature is reduced. This is due to the fact 
that, when scattered from impurities back to their starting point, electron waves 
interfere constructively with their time reversed paths. While this effect is weak at 
high temperatures due to inelastic scattering events, quantum interference becomes 
increasingly important as the temperature is reduced and leads to localization of 
the electrons, albeit on a large length scale; this is generally referred to as "weak 
localization". Indeed, thin metallic films and two-dimensional electron systems 
fabricated on semiconductor surfaces were found to display the predicted logarithmic 
increase of resistivity (Dolan and Osheroff 1979; Bishop et al 1980, 1982; Uren et al 
1980), providing support for the weak localization theory. 

The scaling theory does not explicitly consider the effect of the Coulomb 
interaction between electrons. The strength of the interactions is usually characterized 
by the dimensionless Wigner-Seitz radius, r s = l/(7rn s ) 1 / 2 a_B (here n s is the electron 
density and as is the Bohr radius in a semiconductor). In the experiments mentioned 
above, the Coulomb interactions are relatively weak. Indeed, these experiments are in 
agreement with theoretical predictions (Altshuler, Aronov and Lee 1980) that weak 
electron-electron interactions (r s <C 1) increase the localization even further. As 
the density of electrons is reduced, however, the Wigner-Seitz radius grows and the 
interactions provide the dominant energy of the system. No analytical theory has 
been developed to date in the strongly- interacting limit (r s 3> 1). Finkelstein (1983, 
1984) and Castcllani et al (1984) predicted that for weak disorder and sufficiently 
strong interactions, a 2D system should scale towards a conducting state as the 
temperature is lowered. However, the scaling procedure leads to an increase in the 
effective strength of the interactions and to a divergent spin susceptibility, so that 
the perturbative approach breaks down as the temperature is reduced toward zero. 
Therefore, the possibility of a 2D metallic ground state stabilized by strong electron- 
electron interactions was not seriously considered. 

Recent progress in semiconductor technology has enabled the fabrication of 
high quality 2D samples with very low randomness in which measurements can be 
made at very low carrier densities. The strongly-interacting regime (r s 3> 1) has 
thus become experimentally accessible. Experiments on low-disordered 2D silicon 
samples (Kravchenko et al 1994, 1995, 1996) demonstrated that there are surprising 
and dramatic differences between the behaviour of strongly interacting systems at 
r s > 10 as compared with weakly-interacting systems: with increasing electron 
density, one can cross from the regime where the resistance diverges with decreasing 
temperature (insulating behaviour) to a regime where the resistance decreases strongly 
with decreasing T (metallic behaviour) . These results were met with great scepticism 
and largely overlooked until 1997, when they were confirmed in silicon MOSFETs from 
a different source (Popovic et al 1997) and in other strongly-interacting 2D systems 
(Coleridge et al 1997; Hanein et al 1998a; Papadakis and Shayegan 1998). Moreover, it 
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was found (Simonian et al 1997b; Pudalov et al 1997; Simmons et al 1998) that in the 
strongly-interacting regime, an external magnetic field strong enough to polarize the 
electrons' spins induces a giant positive in-plane magnetoresistance^: and completely 
suppresses the metallic behaviour, thus implying that the spin state is central to the 
high conductance of the metallic state. This finding was in qualitative agreement with 
the prediction of Finkelstcin and Castellani et al that for spin-polarized electrons, only 
an insulating ground state is possible in a disordered 2D system even in the presence 
of strong interactions. Subsequent experiments (Okamoto et al 1999; Kravchenko et 
al 2000b; Shashkin et al 2001, 2002; Vitkalov et al 2001b, 2002; Pudalov et al 
2002b; Zhu et al 2003) have shown that there is a sharp enhancement of the spin 
susceptibility as the metal-insulator transition is approached; indications exist that 
in silicon MOSFETs, the spin susceptibility may actually diverge at some sample- 
independent electron density n x « 8 • 10 10 cm~ 2 . 

In silicon samples with very low disorder potential, the critical density for the 
metal-insulator transition was found to be at or very near n x (Shashkin et al 2001a, 
2002; Vitkalov et al 2001b, 2002), indicating that the metal-insulator transition 
observed in these samples is a property of a clean disorder-free 2D system, rather 
than being a disorder-driven transition. In such samples, the metallic and insulating 
regimes are divided by a temperature-independent separatrix with p s=a 3/i/e 2 , in 
the vicinity of which the resistivity displays virtually universal critical behaviour. 
However, in samples with relatively strong disorder, the electrons become localized at 
densities significantly higher than n x : from 1.44 • 10 11 to 6.6 • 10 11 cm~ 2 (Prus et al 
2002), and even as high as 1.6 • 10 12 cm~ 2 (Pudalov et al 1999), indicating that the 
localization transition in these samples is driven by disorder. 

We suggest that there has been a great deal of confusion and controversy caused by 
the fact that often no distinction has been made between results obtained in systems 
with relatively high disorder and those obtained for very clean samples, and also 
because in many experimental studies, a change in the sign of the derivative di?/dT 
has always been assumed to signal a metal-insulator transition. In this review, we 
focus our attention on results for very clean samples. 

The experimental findings described above were quite unexpected. Once 
accepted, they elicited strong and widespread interest among theorists, with 
proposed explanations that included unusual superconductivity (Phillips et al 1998), 
charging/discharging of contaminations in the oxide (Altshuler and Maslov 1999), 
the formation of a disordered Wigner solid (Chakravarty et al 1999), inter-subband 
scattering (Yaish et al 2000) and many more (for a review, see Abrahams, Kravchenko 
and Sarachik 2001). It is now well-documented that the metallic behaviour in zero 
magnetic field is caused by the delocalizing effect associated with strong electron- 
electron interactions which overpower quantum localization. In the "ballistic regime" 
deep in the metallic state, the conductivity is linear with temperature and derives from 
coherent scattering of electrons by Friedel oscillations (Zala, Narozhny and Aleiner 
2001). Closer to the transition, in the "diffusive" regime, the temperature dependence 
of the resistance is well described by a renormalization group analysis of the interplay 
of strong interactions and disorder (Punnoose and Finkelstcin 2002). Within both 
theories (which consider essentially two limits of the same physical process) an external 
magnetic field quenches the delocalizing effect of interactions by aligning the spins, 



X The fact that the parallel magnetic field promotes insulating behaviour in strongly interacting 2D 
systems was first noticed by Dolgopolov et al (1992). 
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and causes a giant positive magnetoresistance. However, the metal-insulator transition 
itself, as well as the dramatic increase of the spin susceptibility and effective mass in 
its close vicinity, still lack adequate theoretical description; in this region the system 
appears to behave well beyond a weakly interacting Fermi liquid. 

In the next three sections, we describe the main experimental results that 
demonstrate the unexpected presence of a metallic phase in 2D and present an overview 
of recent experiments with emphasis on the anomalous magnetic properties observed 
in the vicinity of the metal-insulator transition. 

2. EXPERIMENTAL RESULTS IN ZERO MAGNETIC FIELD 

2.1. Resistance in zero magnetic field, experimental scaling, reflection symmetry 

The first experiments that demonstrated the unusual temperature dependence of the 
resistivity (Kravchenko et al 1994, 1995, 1996) were performed on low-disordered 
silicon metal-oxide-semiconductor field-effect transistors (MOSFETs) with maximum 
electron mobilities reaching more than 4-10 4 cm 2 /Vs, mobilities that were considerably 
higher than in samples used in earlier investigations. It was the very high quality of the 
samples that allowed access to the physics at electron densities below 10 11 cm" 2 . At 
these low densities, the Coulomb energy, Eq, is the dominant parameter. Estimates 
for Si MOSFETs at n s = 10 11 cm" 2 yield E c « 10 meV, while the Fermi energy, E F , 
is about 0.6 meV (a valley degeneracy of two is taken into account when calculating 
the Fermi energy, and the effective mass is assumed to be equal to the band mass, 
mj.) The ratio between the Coulomb and Fermi energies, r* = Ec/Ep, thus assumes 
values above 10 in these samples. 

Figure ^ (a) shows the temperature dependence of the resistivity measured in 
units of h/e 2 of a high- mobility MOSFET for 30 different electron densities n s varying 
from 7.12 • 10 10 to 13.7 • 10 10 cm~ 2 . If the resistivity at high temperatures exceeds 
the quantum resistance h/e 2 (the curves above the dashed red line), p(T) increases 
monotonically as the temperature decreases, behaviour that is characteristic of an 
insulator. However, for n s above a certain "critical" value, n c (the curves below 
the "critical" curve that extrapolates to 3h/e 2 denoted in red), the temperature 
dependence of p(T) becomes non-monotonic: with decreasing temperature, the 
resistivity first increases (at T > 2 K) and then starts to decrease. At yet higher 
density n s , the resistivity is almost constant at T > 4 K but drops by an order of 
magnitude at lower temperatures, showing strongly metallic behaviour as T — > 0. 

A striking feature of the p(T) curves shown in Figure^^) f° r different n s is that 
they can be made to overlap by applying a (density-dependent) scale factor along the 
T-axis. Thus, the resistivity can be expressed as a function of T/Tq with To depending 
only on n s . This was demonstrated for several samples over a rather wide range of 
electron densities (typically (n c — 2.5 • 10 10 ) < n s < (n c + 2.5 • 10 10 cm -2 )) and in the 
temperature interval 0.2 to 2 K. The results of this scaling are shown in Fig. ^ (b), 
where p is plotted as a function of T/Tq. One can see that the data collapse onto 
two separate branches, the upper one for the insulating side of the transition and the 
lower one for the metallic side. The thickness of the lines is largely governed by the 
noise within a given data set, attesting to the high quality of the scaling. 

The procedure used to bring about the collapse and determine To for each n s was 
the following. No power law dependence was assumed a priori for To versus (n s — n c ); 
instead, the p(T) curves were successively scaled along the T-axis to coincide: the 
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Figure 1. (a): temperature dependence of the B = resistivity in a dilute 
low-disordered Si MOSFET for 30 different electron densities n s ranging from 
7.12 to 13.7- 10 10 cm" 2 , (b): resistivity versus T/T , with T (n s ) chosen to yield 
scaling with temperature. The inset shows the scaling parameter, Tg, versus the 
deviation from the critical point, \n s — n c \; data are shown for silicon MOSFETs 
obtained from three different wafers. Open symbols correspond to the insulating 
side and closed symbols to the metallic side of the transition. From Kravchenko et 
al (1995). 



second "insulating" curve from the top was scaled along the T-axis to coincide with 
the top-most curve and the corresponding scaling factor was recorded, then the third, 
and so on, yielding the upper curve in Fig.^(b) (designated by open symbols). The 
same procedure was applied on the metallic side of the transition starting with the 
highest-density curve, giving the lower curve in Fig. ^(b) (shown as closed symbols). 
A quantitative value was assigned to the scaling factor to obtain To for the insulating 
curves by using the fact that the resistivity of the most insulating (lowest n s ) curve 
was shown (Mason et al , 1995) to exhibit the temperature dependence characteristic 

1 /2~ 

of hopping in the presence of a Coulomb gap: p = po exp (To/T) (Efros and 

Shklovskii 1975). To was determined on the metallic side using the symmetry between 
metallic and insulating curves, as described in more detail in Kravchenko et al (1995). 
For all samples studied, this scaling procedure yields a power law dependence of To 
on \S n \ = |n s — n c \ on both sides of the transition: To oc \S n \^ with the average 
power (3 — 1.60 ± 0.1 for the insulating side and 1.62 ± 0.1 for the metallic side of 
the transition; this common power law can be clearly seen in the inset to Fig.n](b), 
where for each sample the open (insulating side) and filled (metallic side) symbols 
form a single line. The same power law was later observed by Popovic (1997) in 
silicon samples from another source, thus establishing its universality and supporting 
the validity of the scaling analysis. 

Simonian et al (1997a) noted that the metallic and insulating curves are reflection- 
symmetric in the temperature range between approximately 300 mK and 2 K. In 
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Figure 2. (a) For a silicon MOSFET, the normalized resistivity, p* , and 
normalized conductivity, a*, as a function of the gate voltage, V g , at T = 0.35 K. 
Note the symmetry about the line n a = n c . The electron density is given 
by n s = {V g — 0.58V) • 1.1 ■ 10 11 cm -2 , (b) To demonstrate this symmetry 
explicitly, p*(<5 n ) (closed symbols) and cr*(— <5 n ) (open symbols) are plotted as a 
function of 8 n = (n s —n c )/n c . Inset: p*(8 n ) (closed symbols) and cr*(—5 n ) (open 
symbols) versus 8 n at T = 0.3 K and T = 0.9 K, the lowest and highest measured 
temperatures. From Simonian et al (1997a). 



Fig. |21 (a), the normalized resistivity, p* = p/p(n c ), is shown as a function of the 
gate voltage, V g , together with the normalized conductivity, a* = 1/p*; the apparent 
symmetry about the vertical line corresponding to the critical electron density can be 
clearly seen. Fig. [21 (b) demonstrates that the curves can be mapped onto each other 
by reflection, i.e., p*(5 n ) is virtually identical to a*(—6 n ). This mapping holds over a 
range of temperature from 0.3 K to 0.9 K; however, the range |<$ n | over which it holds 
decreases as the temperature is decreased: for example, at T = 0.9 K, p* and a* are 
symmetric for \6 n \ < 0.1, while at T = 0.3 K, they are symmetric only for \8 n \ < 0.05 
(see inset to Fig. 121(b)). Similar symmetry was later reported by Popovic et al (1997) 
and Simmons et al (1998). This implies that there is a simple relation between the 
mechanism for conduction on opposite sides in the vicinity of the transition; the data 
bear a strong resemblance to the behaviour found by Shahar et al (1996) for the 
resistivity near the transition between the quantum Hall liquid and insulator, where it 
has been attributed to charge-flux duality in the composite boson description. It was 
argued by Dobrosavljevic et al (1997) that both the scaling and reflection symmetry 
are consequences of a simple analysis assuming that a T = quantum critical point 
describes the metal-insulator transition. 

2.2. How universal is p(T) ? 

The temperature dependence of p(T) is very similar for clean silicon MOSFETs. As 
an example, Fig. |3 shows the resistivity as a function of temperature of three low- 
disordered samples obtained from different sources. The behaviour is quantitatively 
similar in the critical region in the vicinity of the "separatrix" , which is the horizontal 
curve for which the resistivity is independent of temperature. In all samples, the 
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Figure 3. Universal behaviour in ultra-clean silicon MOSFETs: the resistivity 
versus temperature in three samples from different sources. Note that the 
resistivities are essentially the same at the separatrices for all samples even 
though their critical densities are different, (a) high- mobility sample provided by 
V M Pudalov (graph adopted from Sarachik and Kravchenko 1999), (b) sample 
fabricated by Heemskerk and Klapwijk (1998) (adopted from Kravchenko and 
Klapwijk 2000a) and (c) sample of Heemskerk and Klapwijk but from a different 
wafer (graph adopted from Jaroszyfiski et al 2002). Electron densities in (c) vary 
from 8.55 to 14 ■ 10 10 cm~ 2 (top to bottom). 
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Figure 4. Resistivity versus temperature of a very disordered silicon MOSFET. 
Note that the vertical scale is similar to Fig.|J] The electron densities are (in units 
of 10 11 cm" 2 ): 3.85, 4.13, 4.83, 5.53, 6.23, 7.63, 9.03, 10.4, 11.8, 13.2, 16.0, 18.8, 
21.6, 24.4, 30.0 and 37. Adopted from Pudalov et al (2001). Even though there is 
an apparent crossing point on the p(n a ) isotherms (see the inset), the temperature 
dependence of the resistivity does not resemble the critical behaviour seen in low- 
disordered samples. 

separatrix is remarkably flat at temperatures below 1 K§, and the resistivity is 
essentially the same numerically at slightly less than 3/i/e 2 . The curves below the 
separatrix exhibit strongly metallic temperature dependence (dp /AT > 0) with no 
low-T saturation down to the lowest temperature (40 mK or lower; see Fig|5](b)); the 
drop in resistivity reaches as much as a factor of 10 for the bottom curve in Fig[3] (a). 
Alternative methods used to determine the critical electron density in low-disordered 
samples yield the same value as the density for the separatrix (see sec. 12.3(1 . It is 
important to note that the separatrix in ultra-clean samples represents the "upper 
limit" of the resistivity for which metallic behaviour (as characterized by dp/dT > 0) 
can exist: metallic p(T) has never been observed in any 2D samples at resistivities 
above ~ 3/i/e 2 , in quantitative agreement with the predictions of the renormalization 
group theory (see sec !5.1|) . 

In more disordered samples, however, the behaviour of the resistivity is very 
different. Even though the p(n s ) isotherms apparently cross at some electron density 
(see the inset to FigQJ, the temperature dependence of the resistivity does not resemble 
the critical behaviour seen in low-disordered samples. An example of p(T) curves in 
a disordered sample is shown in Fig0] The sample is insulating at electron densities 
up to ~ 8 ■ 10 11 cm~ 2 which is an order of magnitude higher than the critical density 
in low-disordered samples. The metallic temperature dependence of the resistance 
visible at higher n s does not exceed a few percent (compared to a factor of 10 in 
low-disordered samples). Most importantly, the density corresponding to the crossing 
point shown in the inset to Fig^ldoes not coincide with the critical density determined 
by other methods discussed in the next subsection. 

A metal-insulator transition similar to that seen in clean silicon MOSFETs 
has also been observed in other low-disordered, strongly-interacting 2D systems: p- 

§ At T > 2 K, the resistivity of the separatrix slowly decreases with increasing temperature, as can 
be seen in Fig. Q] where p(T) curves are shown in a much wider temperature interval. 
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Figure 5. For low-disordered 2D hole systems in p-GaAs/AlGaAs, the resistivity 
per square is shown as a function of temperature for B = at various fixed 
hole densities, p. Left hand panel: ISIS (inverted semiconductor-insulator- 
semiconductor) structure with hole densities (from top to bottom) p = 0.89, 0.94, 
0.99, 1.09, 1.19, 1.25, 1.30, 1.50, 1.70, 1.90, 2.50, 3.20, 3.80, 4.50, 5.10, 5.70 and 
6.40-10 10 cm -2 . The inset shows a schematic diagram of the ISIS structure: The 
carriers are accumulated in an undoped GaAs layer situated on top of an undoped 
AlAs barrier, grown over a p+ conducting layer which serves as a back-gate; the 
hole density, p, is varied by applying a voltage to the back gate. From Hanein et 
al (1998a). Right hand panel: Temperature dependence of p in an ultra high 
mobility p-type GaAs/AlGaAs heterostructure at p = 0.48, 0.55, 0.64, 0.72, 0.90, 
1.02, 1.27, 1.98, 2.72 and 3.7210 10 cm" 2 (from top to bottom). From Yoon et al 
(1999). 

type SiGe heterostructures (Coleridge et al 1997), GaAs/AlGaAs heterostructures 
(Hanein et al 1998a; Yoon et al 1999; Mills et al 1999; Noh et al 2002 and others) 
and AlAs heterostructures (Papadakis and Shayegan 1998). It is difficult to make a 
direct comparison of the resistivity observed in different material systems because the 
temperature scales are different, since the Coulomb and Fermi energies depend on the 
effective mass and carrier density. For example, the characteristic temperature below 
which the metallic decrease in the resistivity occurs in p-type GaAs/AlGaAs samples 
is about ten times smaller than in silicon MOSFETs. On the other hand, the values 
of the resistivity are quite similar in the two systems. In Fig. the resistivity is 
plotted as a function of temperature for two p-type GaAs/AlGaAs samples produced 
using different technologies. The data shown in the left hand panel were obtained by 
Hanein et al (1998a) for an inverted semiconductor-insulator-semiconductor (ISIS) 
structure with maximum mobility of /i max = 1-5 • 10 5 cm 2 /Vs, while the right- 
hand panel shows p{T) measured by Yoon et al (1999) on a p-type GaAs/AlGaAs 
heterostructure with peak mobility by a factor of five higher (7 • 10° cm 2 /Vs). The 
interaction parameter, r s , changes between approximately 12 and 32 for the left hand 
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Figure 6. The resistivity as a function of temperature for a disordered p-type 
GaAs/AlGaAs heterostructure at hole densities p = 3.2 — 5.6 • 10 10 cm -2 . From 
Simmons et al (2000). 



plot and from 16 to 44 for the right hand plot 1 1 . In spite of the difference in the 
sample quality and range of densities, the dependence of p(T) on temperature is 
almost the same for the two samples. The main features are very similar to those 
found in silicon MOSFETs: when the resistivity at "high" temperatures exceeds 
the quantum resistance, h/e 2 (i.e., at hole densities below some critical value, p c ), 
the p(T) curves are insulating-like in the entire temperature range; for densities just 
above p c , the resistivity shows insulating-like behaviour at higher temperatures and 
then drops by a factor of 2 to 3 at temperatures below a few hundred mK; and at 
yet higher hole densities, the resistivity is metallic in the entire temperature range. 
Note that the curves that separate metallic and insulating behaviour have resistivities 
that increase with decreasing temperature at the higher temperatures shown; this is 
quite similar to the behaviour of the separatrix in silicon MOSFETs when viewed over 
a broad temperature range (see Fig. (T](a)). However, below approximately 150 mK, 
the separatrix in p-type GaAs/AlGaAs heterostructures is independent of temperature 
(Hanein et al , 1998b), as it is in Si MOSFETs below approximately 2 K. The resistivity 
of the separatrix in both systems extrapolates to 2 — 3 h/e 2 as T — > 0, even though 
the corresponding carrier densities are very different. 

As in the case of highly disordered silicon MOSFETs, no critical behaviour of 
resistance is observed in disordered GaAs/AlGaAs heterostructures. An example is 
shown in Fig. where the temperature dependence of the resistivity at B = is 
plotted for hole densities p = 3.2 — 5.6T0 10 cm~ 2 . Monotonic localized behaviour is 
observed even when the "high-temperature" resistivity lies well below h/e 2 , at carrier 
densities up to 4.6T0 10 cm -2 ; both samples shown in the previous figure would be in 
the deeply metallic regime at this hole density. Above this density, the decrease in 
resistivity with decreasing temperature is very small (about 10% in the temperature 
interval 0.7 to 0.1 K). 

|| These r 3 values were calculated assuming that the effective mass is independent of density and 
equal to 0.37 m e , where m e is the free-electron mass. 
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As indicated by the experimental results presented above, the p(T) curves are 
nearly universal in the vicinity of the metal-insulator transition, but only in samples 
with very weak disorder potential. The strength of the disorder is usually characterized 
by the maximum carrier mobility, p, max . In general, the higher the maximum mobility 
(i.e. the lower the disorder), the lower the carrier density at which the localization 
transition occurs. This was shown empirically by Sarachik (2002) to hold over a 
broad range (five decades in density) for all materials studied: the critical density 
follows a power law dependence on peak mobility (or scattering rate). However, the 
data exhibit some scatter, and the correlation is not exact. Thus, for example, the 
peak hole mobilities in samples used by Hanein et al (1998a) and Simmons et al 
(1998) are similar, while the localization transition occurs in the latter at a value of p 
several times higher than the former. This may be due to sample imperfections (e.g., 
a slightly inhomogeneous density distribution), which are important at low carrier 
densities, while the maximum mobility is reached at relatively high carrier densities 
and may therefore be relatively insensitive to such effects. 

A better indicator of the strength of the disorder potential near the MIT is 
how low the carrier density is at which the localization transition occurs. In silicon 
MOSFETs, the experimental results obtained to date suggest that the resistivity near 
the transition approaches universal behaviour for samples in which the transition to 
a strongly localized state occurs at n s < 1 ■ 10 11 cm~ 2 . (In p-type GaAs/AlGaAs 
heterostructures, the corresponding density separating universal and non-universal 
behaviour appears to be about an order of magnitude lower, although the data are 
currently insufficient to determine the value reliably.) Below, we will argue that 
the "universal" metal-insulator transition in very clean samples is not driven by 
disorder but by some other mechanism, possibly of magnetic origin. In contrast, 
the transition is not universal in more disordered samples and is presumably due to 
Anderson localization, which is strong enough to overpower the metallic behaviour at 
low densities. 

2.3. Critical density 

To verify whether or not the separatrix corresponds to the critical density, an 
independent determination of the critical point is necessary: comparison of values 
obtained using different criteria provides an experimental test of whether or not 
a true MIT exists at B = 0. One obvious criterion, hereafter referred to as 
the "derivative criterion" , is a change in sign of the temperature derivative of the 
resistivity, dp/dT; this is the criterion often used to identify the MIT. A positive 
(negative) sign of the derivative at the lowest achievable temperatures is empirically 
associated with a metallic (insulating) phase. A weakness of this criterion is that it 
requires extrapolation to zero temperature. A second criterion can be applied based 
on an analysis of a temperature-independent characteristic, namely, the localisation 
length L extrapolated from the insulating phase. These two methods have been applied 
to low-disordered silicon MOSFETs by Shashkin et al (2001b) and Jaroszyhski et al 
(2002). 

As mentioned earlier, the temperature dependence of the resistance deep in the 
insulating phase obeys the Efros-Shklovskii variable-range hopping form (Mason et al 
1995); on the other hand, closer to the critical electron density at temperatures that 
are not too low, the resistance has an activated form p cx e Ea ^ kBT (Pepper et al 1974; 
Pudalov et al 1993; Shashkin et al 1994) due to thermal activation to the mobility 
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Figure 7. Activation energy (diamonds) and square root of the threshold voltage 
(circles) versus electron density in zero magnetic field in a low-disordered silicon 
MOSFET. The inset shows current-voltage characteristics recorded at pa 30 and 
210 mK, as labelled; note that the threshold voltage is essentially independent of 
temperature. From Shashkin et al (2001b). 
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Figure 8. Conductivity versus temperature for different n s in a highly disordered 
silicon MOSFET. (The error bars show the size of the fluctuations of the 
conductivity with time.) The critical electron densities obtained by the derivative 
and localization length criteria (n* and n c i, correspondingly) are marked by 
arrows. n g is the density corresponding to the onset of glassy behaviour; see 
text. The figure is adopted from Bogdanovich and Popovic (2002). 
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edge. Fig. [7| shows the activation energy E a as a function of the electron density 
(diamonds); the data can be approximated by a linear function which yields, within 
the experimental uncertainty, the same critical electron density as the "derivative 
criterion" . 

The critical density can also be determined by studying the nonlinear current- 
voltage I — V characteristics on the insulating side of the transition. A typical low- 
temperature I — V curve is close to a step-like function: the voltage rises abruptly at 
low current and then saturates, as shown in the inset to Fig. the magnitude of the 
step is 2V C . The curve becomes less sharp at higher temperatures, yet the threshold 
voltage, V c , remains essentially unchanged. Closer to the MIT, the threshold voltage 
decreases, and at n s = n c \ = 0.795 • 10 11 cm" 2 , the I — V curve is strictly linear 
(Shashkin et al 2001b). According to Polyakov and Shklovskii (1993) and Shashkin et 
al (1994), the breakdown of the localized phase occurs when the localized electrons at 
the Fermi level gain enough energy to reach the mobility edge in an electric field, V c /d, 
over a distance given by the localization length, L, which is temperature- independent: 

eV c (n s ) L(n s )/d = E a (n s ) 

1/2 

(here d is the distance between the potential probes). The dependence of V c (n s ) 
near the MIT is linear, as shown in Fig.0by closed circles, and its extrapolation to zero 
threshold value again yields approximately the same critical electron density as the 

1/2 

two previous criteria. The linear dependence V c 1 (n s ), accompanied by linear E a (n s ), 
signals the localization length diverging near the critical density: L(n s ) oc 1/ (n c — n s ). 

These experiments indicate that in low-disordered samples, the two methods — 
one based on extrapolation of p(T) to zero temperature and a second based on the 
behaviour of the temperature-independent localization length — give the same critical 
electron density: n c ss n c \. This implies that the separatrix remains "flat" (or 
extrapolates to a finite resistivity) at zero temperature. Since one of the methods is 
independent of temperature, this equivalence supports the existence of a true T = 
MIT in low-disordered samples in zero magnetic field. 

In contrast, we note that in highly-disordered samples, the localization length 
method yields a critical density noticeably lower than the derivative criterion. Figure|H] 
shows the (time-averaged) conductivity (a) as a function of T for different n s . d(a) / dT 
changes sign at electron density n* = 12.9 • 10 11 cm" 2 . The activation energy, 
however, vanishes at n c \ rs 5.2 • 10 11 cm" 2 , which is more than a factor of two lower 
than n* s : at densities between these two values, the resistivity does not diverge as 
T — > 0, even though it exhibits an insulating-like temperature dependence. Thus, 
these two different methods yield different "critical densities" for a sample with strong 
disorder. Moreover, from the study of low-frequency resistance noise in dilute silicon 
MOSFETs, Bogdanovich and Popovic (2002) and Jaroszyhski et al (2002) have found 
that the behaviour of several spectral characteristics indicates a dramatic slowing 
down of the electron dynamics at a well-defined electron density n g , which they have 
interpreted as an indication of a (glassy) freezing of the 2D electron system. In low- 
disordered samples, n g nearly coincides with n c , while in highly-disordered sample, 
n g lies somewhere between n c \ and n*, as indicated in Fig. |H1 The width of the 
glass phase (n c i < n s < n g ) thus strongly depends on disorder, becoming extremely 
narrow (or perhaps even vanishing) in low-disordered samples. The strong dependence 
on disorder of the width of the metallic glass phase is consistent with predictions of 
the model of interacting electrons near a disorder-driven metal-insulator transition 
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Figure 9. (a-d) The left hand panels show data for the resistivity versus 
temperature at B = in a disordered p-type GaAs/AlGaAs quantum well, 
illustrating the transition from insulating to metallic behaviour as the density 
increases. The right hand panels show the corresponding magnetoresistance traces 
for temperatures of 147, 200, 303, 510, 705, and 910 mK. From Simmons et al 
(2000). 

(Pastor and Dobrosavljevic 1999). These observations all suggest that the origin of 
the metal-insulator transition is different in clean and strongly-disordered samples. 

2-4- Does weak localization survive in the presence of strong interactions? 

The theory of weak localization was developed for noninteracting systems, and it was 
not a priori clear whether it would work in the presence of strong interactions. In 
2000, Simmons et al studied transport properties of a dilute modulation-doped p- 
type GaAs/AlGaAs quantum well and observed a temperature-dependent negative 
magnetoresistance, consistent with the suppression of the coherent backscattering by 
the perpendicular magnetic field. Magnetoresistance curves obtained by Simmons et 
al are plotted in the right hand panel of Fig. A characteristic peak develops 
in the resistivity at B — as the temperature is decreased, signalling that the 
weak localization is still present at p as low as 4.5 • 10 10 cm~ 2 , corresponding 
to the interaction parameter r s ~ 15. Simmons et al successfully fitted their 
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Figure 10. The left hand panels show resistance per square as a function 
of temperature for 2D holes in an ultra-clean GaAs/AlGaAs heterostructure for 
various values of the gate bias. The right hand panels show the variation of 
longitudinal resistance with perpendicular magnetic field for the same sample at 
various gate biases. The solid grey lines represent the fit described in the text. 
From Mills et al (2001). 



magnetoresistance data by the Hikami-Larkin formula (Hikami et al 1980) 

and obtained reasonable values of the phase-breaking time, ~ 10 to 30 ps, in 
the temperature interval 1 to 0.15 K (here r is the elastic scattering time, is 
the Digamma function, tb — h/4eB±D, and D is the diffusion coefficient). Weak 
negative magnetoresistance was also observed by Brunthaler et al (2001) in silicon 
MOSFETs at electron densities down to 1.5 • 10 11 cm -2 , although they found values 
of T0 that were about an order of magnitude shorter than those expected theoretically, 
~ ah 2 /e 2 kBT. 

However, the agreement with non-interacting theory breaks down at higher 
interaction strengths. Mills et al (2001) studied p- type GaAs/AlGaAs heterostructures 
of much higher quality which remained metallic down to p w 3 • 10 9 cm~ 2 
(corresponding to r s ~ 60 provided the effective mass does not change), and found 
the coherent backscattering to be almost completely suppressed at these ultra-low hole 
densities. In the right hand panel of Fig. the magnetoresistance traces are shown 
at T w9 mK for four different carrier densities. The width of the characteristic peak 
at B± = 0, visible in the top two curves, is approximately as expected from the theory, 



Metal-insulator transition in 2D 



17 



10 



3, 
b 




-0.4 



-0.2 



0.2 

(tesla) 



0.90 
0.87 
0.82 

0.4 



-(b) 


e 2 /h 123^ 







X 1 

x6 
x50 



-0.2 



-0.1 



0.1 

(tesla) 



0.2 



Figure 11. Longitudinal magnetoconductivity of a low-disordered silicon 
MOSFET in a weak perpendicular magnetic field at T = 42 mK for a range 
of electron densities indicated near each curve in units of 10 11 cm -2 . The curves 
in (b) are vertically shifted and the two lowest curves are multiplied by 6 (the 
middle one) and 50 (the lower one). The thick dashed line in (b) is a fit by Eg. HI 
with b = 0.6 and = 30 ps. From Rahimi et al (2003). 



but its magnitude is about a factor of 30 smaller than expected. At slightly higher p 
(the two bottom curves), the peak is not seen at all. 

In principle, the strong disagreement between the expected and measured peak 
magnitudes might be due to the fact that the theory for coherent backscattering is not 
applicable to a system with resistivity of the order of h/e 2 . However, this is not the 
source of the disagreement, as the resistivity in the experiments of Mills et al (2001) 
is in the same range as in the experiments of Simmons et al (2000) (cf Figs. I§1 and lTU|) . 
Therefore the suppression of the peak is apparently related to stronger interactions 
(higher r s ) in the Mills et al samples rather than to high values of p. 

The left hand panel of Fig. ^| shows the temperature dependence of the B = Q 
resistance in the ultra low density sample of Mills et al over the temperature range 5 
to 100 mK. The dashed lines show the expected corrections to the resistivity caused by 
weak localization calculated using the equation Act(T) = b e 2 /h ln^/Y), where b is 
a constant expected to be universal and equal to 1/tt. The calculated dependence 
is clearly at variance with the measurements; rather, at very low temperatures, 
the resistance becomes nearly constant. Fitting the theoretical expression to the 
experimental data, Mills et al found that the upper limits for b are from one to 
nearly two orders of magnitude smaller than the b = 1/tt expected from the theory. 

The disappearance of weak localization corrections near the MIT has also been 
observed by Rahimi et al (2003) in low-disordered silicon MOSFETs. The results are 
shown in Fig. 1111 At higher n s (the upper curves in Fig Jlll fa) ), the characteristic dip 
is observed in the magnetoconductance at zero magnetic field. As follows from Eq. ^ 
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the magnitude of the dip is expected to be equal to (be 2 g v /h) ln^/r), and should 
therefore exhibit a weak (double-logarithmic) increase as the average conductivity 
decreases provided the variations in electron density are small, as they are in this case. 
This is not what is observed in the experiment: as one approaches the transition, the 
magnitude of the dip decreases sharply, and at the critical electron density (the lowest 
curve in Fig. (a)), the dip is no longer seen on the scale of this figure. However, 
the shape of the magnetoconductivity does not change significantly with decreasing 
n s as illustrated by the middle curve in Fig. 1111 (b). which shows <j(B±) multiplied by 
six (n s — 0.85 • 10 11 cm -2 ) to make it quantitatively similar to the upper curve. This 
similarity demonstrates that the functional form of the a(B±) dependence, described 
by the expression in square brackets in Eq. does not change noticeably as the 
density is reduced from 1.23 • 10 11 to 0.85 • 10 11 cm~ 2 ; instead, it is the magnitude 
of the effect that rapidly decreases upon approaching the MIT. At yet lower density, 
n s = 0.82 ■ 10 11 cm' 2 , the magnitude of the dip does not exceed 2% of that for 
n s = 1.23 • 10 11 cm~ 2 (compare the upper and the lower curves in Fig. 1111 (b)). 

It may seem surprising that a change in n s by only a factor of 1.5 (from 
n s — 1.23 ■ 10 11 to n s = 0.82 • 10 11 cm -2 ) results in such a dramatic suppression of the 
quantum localization. It is interesting to note in this connection that the experimental 
data on silicon MOSFETs described in sections 14.21 and [4.31 reveal a sharp increase of 
the effective mass in the same region of electron densities where the suppression of 
the weak localization is observed. Due to the strong renormalization of the effective 
mass, the ratio between the Coulomb and Fermi energies, r* = 2 r s (m*/mb), grows 

— 1/2 

much faster than n s reaching values greater than 50 near the critical density. 

The apparent absence of localization at and just above the critical density may 
account for the existence of a flat separatrix at n s = n c (see Fig. [3J ■ If the localization 
were present, the temperature- independent curve would require that the temperature 
dependence of p due to localization be cancelled exactly over a wide temperature range 
by a temperature dependence of opposite sign due to interactions, a coincidence which 
seems very improbable for two unrelated mechanisms. Note that at resistivity levels 
of order or greater than h/e 2 , the quantum corrections to the resistivity are expected 
to be very strong and cannot be easily overlooked. The calculated temperature 
dependence of the resistivity, expected for non-interacting electrons (Abrahams et al 
1979), is shown in FigEl by the dashed curve: at 100 mK, quantum localization 
is expected to cause a factor of more than 30 increase in resistivity, in strong 
contradiction with the experiment which shows it to be constant within ±5%. 

3. THE EFFECT OF A MAGNETIC FIELD 

3.1. Resistance in a parallel magnetic field 

In ordinary metals, the application of a parallel magnetic field (Bu) does not lead 
to any dramatic changes in the transport properties: if the thickness of the 2D 
electron system is small compared to the magnetic length, the parallel field couples 
largely to the electrons' spins while the orbital effects are suppressed. Only weak 
corrections to the conductivity are expected due to electron-electron interactions (Lee 
and Ramakrishnan 1982, 1985). It therefore came as a surprise when Dolgopolov et 
al (1992) observed a dramatic suppression of the conductivity in dilute Si MOSFETs 
by a parallel in-plane magnetic field B\\ . The magnetoresistance in a parallel field 
was studied in detail by Simonian et al (1997b) and Pudalov et al (1997), also in Si 



Metal-insulator transition in 2D 



19 



1(T 



^> 1 

10 



Q. 



10 L 



f 



y localization 
\ ^ theory 
\ 

\ 

\ 

\ 

experiment 
Nl N W H* * r * 



r/r F = o.oi 



0.0 0.2 0.4 0.6 0.8 1.0 

T(K) 



Figure 12. Resistivity at the separatrix in a low-disordered silicon MOSFET as 
a function of temperature (the solid curve) compared to that calculated from the 

one-parameter scaling theory using d ^ np ( L <t>} = — /3(p) (here La oc p 'T — p / 2 

d in Lfi ' ™ 

is the phase-breaking length, /3(p) is the scaling function approximated by 

f3(p) = — ln(l + ap) following Altshuler et al (2000), a = 2/tt, p is measured 

in units of h/e 2 , and p and 7 are constants equal to 3 and 0.5 respectively). 

As the dashed line shows, the resistivity of a "conventional" (noninteracting) 2D 

system should have increased by a factor of more than 30 when the temperature 

has been reduced to 100 mK. From Kravchcnko and Klapwijk (2000a). 



MOSFETs. In the left hand part of Fig. ^3 the resistivity is shown as a function of 
parallel magnetic field at a fixed temperature of 0.3 K for several electron densities. 
The resistivity increases sharply as the magnetic field is raised, changing by a factor 
of about 4 at the highest density shown and by more than an order of magnitude 
at the lowest density, and then saturates and remains approximately constant up to 
the highest measuring field, Bn — 12 tesla. The magnetic field where the saturation 
occurs, -Bsat, depends on n s , varying from about 2 tesla at the lowest measured density 
to about 9 tesla at the highest. The metallic conductivity is suppressed in a similar 
way by magnetic fields applied at any angle relative to the 2D plane (Kravchenko et al 
1998) independently of the relative directions of the measuring current and magnetic 
field (Simonian et al 1997b; Pudalov et al 2002a). All these observations suggest that 
the giant magnetoresistance is due to coupling of the magnetic field to the electrons' 
spins. Indeed, from an analysis of the positions of Shubnikov-de Haas oscillations in 
tilted magnetic fields, Okamoto et al (1999) and Vitkalov et al (2000, 2001a) have 
concluded that in MOSFETs at relatively high densities, the magnetic field i? sa t is 
equal to that required to fully polarize the electrons' spins. 

In p-type GaAs/AlGaAs heterostructures, the effect of a parallel magnetic field 
is qualitatively similar, as shown in the right hand part of Fig. 1131 The dependence 
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Figure 13. Left hand panel: resistivity versus parallel magnetic field measured 
at T = 0.29 K on low-disordered silicon sample. Different symbols correspond to 
densities from 1.01 to 2.17- 10 11 cm" 2 (adopted from Pudalov et al 1997). Right 
hand panel: resistivity as a function of By in a p-GaAs/AlGaAs heterostructure 
at 50 mK at the following hole densities, from the bottom: 4.11, 3.23, 2.67, 
2.12, 1.63, 1.10, 0.98, 0.89, 0.83, 0.79, 0.75, 0.67-10 10 cm" 2 . The solid lines are 
for hole densities above p c and the open circles are for densities below p c . The 
solid circles denote the experimentally determined critical magnetic fields, and the 
dashed line is a guide to the eye. i?|J, the boundary separating the high and the 
low field regions, is shown by the dotted line. Adopted from Yoon et al (2000). 
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Figure 14. Resistivity versus temperature for five different fixed magnetic fields 
applied parallel to the plane of a low-disordered silicon MOSFET. The electron 
density is 8.83 ■ 10 10 cm" 2 . From Simonian et al (1997b). 
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of p on -By does not saturate to a constant value as in Si MOSFETs, but continues to 
increase with increasing field, albeit at a considerably slower rate. This is presumably 
due to strong coupling of the parallel field to the orbital motion arising from the 
finite layer thickness (see Das Sarma and Hwang 2000), an effect that is more 
important in GaAs/AlGaAs heterostructures than in silicon MOSFETs because of 
a much thicker layer. As in the case of Si MOSFETs, there is a distinct knee that 
serves as a demarcation between the behaviour in low and high fields. For high 
hole densities, Shubnikov-de Haas measurements (Tutuc et al 2001) have shown that 
this knee is associated with full polarization of the spins by the in-plane magnetic 
field. However, unlike Si MOSFETs, the magnetoresistance in p-GaAs/AlGaAs 
heterostructures has been found to depend on the relative directions of the measuring 
current, magnetic field, and crystal orientation (Papadakis et al 2000); one should note 
that the crystal anisotropy of this material introduces added complications. In p-SiGe 
heterostructures, the parallel field was found to induce negligible magnetoresistance 
(Senz et al 1999) because in this system the parallel field cannot couple to the spins 
due to very strong spin-orbit interactions. 

Over and above the very large magnetoresistance induced by an in-plane magnetic 
fields, an even more important effect of a parallel field is that it causes the zero- field 2D 
metal to become an insulator (Simonian et al 1997b; Mertes et al 2001; Shashkin et 
al 2001b; Gao et al 2002). Figure IHI shows how the temperature dependence of 
the resistance changes as the magnetic field is increased. Here, the resistivity of a 
Si MOSFET with fixed density on the metallic side of the transition is plotted as a 
function of temperature in several fixed parallel magnetic fields between and 1.4 tesla. 
The zero-field curve exhibits behaviour typical for "just-metallic" electron densities: 
the resistivity is weakly-insulating at T > T max « 2 K and drops substantially as 
the temperature is decreased below T max . In a parallel magnetic field of only 1.4 tesla 
(the upper curve), the metallic drop of the resistivity is completely suppressed, so that 
the system is now strongly insulating in the entire temperature range. The effect of 
the field is negligible at temperatures above T max , i.e., above the temperature below 
which the metallic behaviour in B = sets in. 

The extreme sensitivity to parallel field is also illustrated in Fig. ^| where the 
temperature dependence of the resistivity is compared in the absence (a) and in the 
presence (b) of a parallel magnetic field. For B\\ = 0, the resistivity displays the 
familiar, nearly symmetric (at temperatures above 0.2 K) critical behaviour about the 
separatrix (the dashed line). However, in a parallel magnetic field of B\\ = 4 tesla, 
which is high enough to cause full spin polarization at this electron density, all the 
p(T) curves display "insulating-like" behaviour, including those which start below 
h/e 2 at high temperatures. There is no temperature-independent separatrix at any 
electron density in a spin-polarized electron system (Simonian et al 1997b; Shashkin et 
al 2001b). 

This qualitative difference in behaviour demonstrates convincingly that the spin- 
polarized and unpolarized states behave very differently. This rules out explanations 
which predict qualitatively similar behaviour of the resistance regardless of the 
degree of spin polarization. In particular, the explanation of the metallic behaviour 
suggested by Das Sarma and Hwang (1999) (see also Lilly et al 2003), based on the 
temperature-dependent screening, predicts metallic-like p(T) for both spin-polarized 
and unpolarized states, which is in disagreement with experiment (for more on this 
discrepancy, see Mertes et al 2001). 
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Figure 15. Temperature dependence of the resistivity of a low-disordered silicon 
MOSFET at different electron densities near the MIT in zero magnetic field (a), 
and in a parallel magnetic field of 4 tesla (b). The electron densities are indicated 
in units of 10 11 cm~ 2 . Dashed curves correspond to n a = n c \ which is equal to 
0.795 • I0 11 cm~ 2 in zero field and to 1.155 • 10 11 cm~ 2 in By = 4 tesla. From 
Shashkin et al (2001b). 



3.2. Scaling of the magnetoresistance; evidence for a phase transition 

There have been many attempts to obtain a quantitative description of the 
magnetoresistance as a function of the carrier density and temperature over the entire 
field range, including the saturation region. Attempts to obtain a collapse of the 
magnetoresistance onto a single scaled curve have yielded scaling at either low or high 
magnetic field; over a wide range of temperatures but only at the metal-insulator 
transition; or in a wide range of carrier densities, but only in the limit of very low 
temperatures. Simonian et al (1998) found that at the transition, the deviation of 
the magnetoconductivity from its zero-field value, Act = cr(Bu) — cr(0), is a universal 
function of B\\/T. Although the quality of the scaling is good, it breaks down rather 
quickly as one moves into the metallic phase. 

Two scaling procedures have been recently proposed; although they differ in 
procedure and yield results that differ somewhat in detail, the major conclusions 
are essentially the same, as described below, and imply that there is an approach to 
a quantum phase transition at a density near n c . 

Shashkin et al (2001a) scaled the magnetoresistivity in the spirit of the theory 
of Dolgopolov and Gold (2000), who predicted that at T = 0, the normalized 
magnetoresistance is a universal function of the degree of spin polarization, P = 
g* /ibB\\/2Ef — g*m*HBB\\/nh 2 n s (here m* is the effective mass and g* is the g- 
factor) . Shashkin et al scaled the data obtained in the limit of very low temperatures 
where the magnetoresistance becomes temperature-independent and, therefore, can be 
considered to be at its T = value. In this regime, the normalized magnetoresistance, 




Figure 16. (a) Low-temperature magnetoresistance of a clean silicon MOSFET 
in parallel magnetic field at different electron densities above n c . (b) Scaled 
curves of the normalized magnetoresistance versus Bn/B a . The electron densities 
are indicated in units of 10 11 cm -2 . Also shown by a thick dashed line is 
the normalized magnetoresistance calculated by Dolgopolov and Gold (2000). 
Adopted from Shashkin et al (2001a). 



p(-B||)/p(0), measured at different electron densities, collapses onto a single curve 
when plotted as a function of B\\/B c (here B c is the scaling parameter, normalized 
to correspond to the magnetic field B sat at which the magnetoresistance saturates). 
An example of how p(B^), plotted in Fig. ^] (a), can be scaled onto a universal 
curve is shown in Fig. 1161 fb). The resulting function is described reasonably well by 
the theoretical dependence predicted by Dolgopolov and Gold. The quality of the 
scaling is remarkably good for B\\fB c < 0.7 in the electron density range 1.08 • 10 11 to 
10 • 10 11 cm -2 , but it breaks down as one approaches the metal-insulator transition 
where the magnetoresistance becomes strongly temperature-dependent even at the 
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Figure 17. Scaling parameter B c (corresponding to the field required for full 
spin polarization) as a function of the electron density. An expanded view of the 
region near n c is displayed in the inset. From Shashkin et al (2001a). 



lowest experimentally achievable temperatures. As shown in Fig. El the scaling 
parameter is proportional over a wide range of electron densities to the deviation 
of the electron density from its critical value: B c oc (n s — n c ). 

Vitkalov et al (2001b) succeeded in obtaining an excellent collapse of 
magnetoconductivity data over a broad range of electron densities and temperatures 
using a single scaling parameter. They separated the conductivity into a field- 
dependent contribution, o~(B\\) — er(oo), and a contribution that is independent 
of magnetic field, a (00). The field-dependent contribution to the conductivity, 
cr(0) — o~(Bu), normalized to its full value, c(0) — a (00), was shown to be a universal 
function of B/B a : 

M = ^,/B.) (2) 
er(0) — a (00) 

where B a (n Sl T) is the scaling parameter. Applied to the magnetoconductance curves 
shown in Fig.^|(a) for different electron densities, the above scaling procedure yields 
the data collapse shown in Fig. ^| (b) . 

Remarkably, similar scaling holds for curves obtained at different temperatures. 
This is demonstrated in Fig. 1191 which shows the scaled magnetoconductance measured 
at a fixed density and different temperatures. Altogether, the scaling holds for 
temperatures up to 1.6 K over a broad range of electron densities up to 4n c . 

Fig. EH shows the scaling parameter B a plotted as a function of temperature 
for different electron densities n s > n c . The scaling parameter becomes smaller as 
the electron density is reduced; for a given density, B a decreases as the temperature 
decreases and approaches a value that is independent of temperature, B a (T = 0). 
As the density is reduced toward n c , the temperature dependence of B a dominates 
over a broader range and becomes stronger, and the low-temperature asymptotic 
value becomes smaller. Note that for electron densities below 1.36 • 10 11 cm" 2 , B a 
is approximately linear with temperature at high T. The behaviour of the scaling 
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Figure 18. (a) Conductivity of a low-disordered silicon sample versus in-plane 
magnetic field at different electron densities in units of 10 11 cm -2 , as labelled; 
T = 100 mK. (b) Data collapse obtained by applying the scaling procedure, Eg 121 
to the curves shown in (a). Adopted from Vitkalov et al (2001b). 

parameter B a {T) can be approximated by an empirical fitting function: 

B a (n s ,T) = A(n s )[A(n s ) 2 + T 2 ] 1/2 - 

The solid lines in Fig. [201 ( a ) are fits to this expression using A(n s ) and A(n s ) as 
fitting parameters. As can be inferred from the slopes of the curves of Fig. [201 (a), 
the parameter A(n s ) is constant over most of the range and then exhibits a small 
increase (less than 20%) at lower densities. As shown in Fig. [201 (b) , the parameter 
A decreases sharply with decreasing density and extrapolates to zero at a density no 
which is within 10% of the critical density n c m 0.85 -10 11 cm" 2 for the metal-insulator 
transition. 

The scaling parameters B c and A in the analysis by Shashkin et al (2001a) and 
Vitkalov et al (2001b) represent energy scales \ibB c and fcsA, respectively, which 
vanish at or near the critical electron density for the metal-insulator transition. At 
high electron densities and low temperatures T < ^bBc/Ub (corresponding to T < A), 
the system is in the zero temperature limit. As one approaches n c , progressively lower 
temperatures are required to reach the zero temperature limit. At n s = hq, the 
energies ^bB c and fc B A vanish; the parameter B a — > as T — > 0; the system thus 
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Figure 19. Data collapse obtained by applying the scaling procedure, Eq|2l 
to the magnetoconductivity at different temperatures for electron density n s = 
9.4 • 10 10 cm -2 . The inset shows the conductivity at different temperatures as a 
function of magnetic field. Adopted from Vitkalov et al (2001b). 



exhibits critical behaviour (Sondhi et al 1997; Vojta 2003), signalling the approach to 
a new phase in the limit T = at a critical density sa n c . 

4. SPIN SUSCEPTIBILITY NEAR THE METAL-INSULATOR 
TRANSITION 

4-1- Experimental measurements of the spin susceptibility 

In the Fermi-liquid theory, the electron effective mass and the g-factor (and, 
therefore, the spin susceptibility x ^ 9* m *) ar e renormalized due to electron-electron 
interactions (Landau 1957). Earlier experiments (Fang and Stiles 1968; Smith and 
Stiles 1972), performed at relatively small values of r s ~ 2 to 5, confirmed the expected 
increase of the spin susceptibility. More recently, Okamoto et al (1999) observed a 
renormalization of \ by a factor of ~ 2.5 at r s up to about 6 (see Fig.|2H (a)). At yet 
lower electron densities, in the vicinity of the metal-insulator transition, Kravchenko et 
al (2000b) have observed a disappearance of the energy gaps at "cyclotron" filling 
factors which they interpreted as evidence for an increase of the spin susceptibility by 
a factor of at least 5. 

It was noted many years ago by Stoner that strong interactions can drive an 
electron system toward a ferromagnetic instability (Stoner 1946). Within some 
theories of strongly interacting 2D systems (Finkelstein 1983, 1984; Castellani et 
al 1984), a tendency towards ferromagnetism is expected to accompany metallic 
behaviour. The experiments discussed in Section 13.21 which indicate that B c and A 
vanish at a finite density prompted Shashkin et al (2001a) and Vitkalov et al (2001b) 
to suggest that spontaneous spin polarization may indeed occur at or near the critical 
electron density in the limit T = 0. The data obtained in these experiments provide 
information from which the spin susceptibility, x, can be calculated in a wide range of 
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Figure 20. (a) B a versus temperature for different electron densities; the solid 
lines are fits to the empirical form B^(n s ,T) = A(n a )[[A(n s )] 2 +T 2 ] 1 / 2 . (b) The 
parameter A as a function of electron density; the solid line is a fit to the critical 
form A = A (n s - no)". From Vitkalov et al (2001b). 

densities. In the clean limit, the band tails are small (Vitkalov et al 2002; Dolgopolov 
and Gold 2002; Gold and Dolgopolov 2002) and can be neglected, and the magnetic 
field required to fully polarize the spins is related to the ^-factor and the effective mass 
by the equation g* /ibB c = 2Ef = irh 2 n s /m* (here, the two-fold valley degeneracy in 
silicon has been taken into account). Therefore, the spin susceptibility, normalized by 
its "non-interacting" value, can be calculated as 

X g*m* nh n s 
Xo gom b 2n B B c m b ' 

Fig. 1211 shows the normalized spin susceptibility (Shashkin et al 2001a) and its 
inverse (Vitkalov et al 2002) obtained using the above expression. The values deduced 
by both groups indicate that g*m* diverges ((<?*m*) _1 extrapolates to zero) in silicon 
MOSFETs at a finite density close or equal to n c . Also shown on the same figure 
are the data of Pudalov et al obtained from an analysis of Shubnikov-de Haas (SdH) 
measurements in crossed magnetic fields (see section E?2*| . The susceptibilities obtained 
by all three groups on different samples, by different methods and in different ranges of 
magnetic field, are remarkably similar (on the mutual consistency of the data obtained 
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Figure 21. Upper graph: normalized spin susceptibility vs n s obtained from 
the data in Fig. 1171 The dashed line is a guide to the eye. The vertical 
dashed line denotes the position of the critical density for the metal-insulator 
transition. Lower graph: the inverse of the normalized spin susceptibility Xo/x* 
versus electron density obtained by Vitkalov et al (2001b), plotted with data of 
Shashkin et al (2001a) and Pudalov et al (2002b). The solid curve is a fit to 
the critical form xo/x* = A(n s — no) a for the data of Vitkalov et al (2001b) 
(excluding the point shown at xo/x* = 0)- Adopted from Vitkalov et al (2002). 



on different samples by different groups, see also Kravchenko et al (2002) and Sarachik 
and Vitkalov (2003)). 

A novel and very promising method has recently been used by Prus et al (2003) 
to obtain direct measurements of the thermodynamic spin susceptibility. The method 
entails modulating the (parallel) magnetic field by an auxiliary coil and measuring the 
AC current induced between the gate and the 2D electron gas, which is proportional 
to dfi/dB (where fj, is the chemical potential). Using Maxwell's relation, 

_ dM 
dB dn s ' 
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Figure 22. Density-dependence of m* g* in an ultra-clean 2D electron system 
in GaAs/AlGaAs determined by two different methods. The solid data points 
are obtained from tilted-field Shubnikov-de Haas measurements with different 
configurations of Landau levels. The parallel dashed lines indicate a power law 
dependence of m* g* with a single exponent for all level configurations. The 
open circles show nonmonotonic density-dependence of m*g* derived from the 
full polarization field, B c , using the parallel field method. The inset shows the net 
spin for n s = 1- 10 10 cm~ 2 with interpolated regime (solid line) and extrapolated 
regime (dotted line). B C =4.9T from the in-plane field method and B ex t=6.5 T 
from extrapolation of the tilted-field method. The thick solid line represents 
extrapolation of m* g* to the P=0 limit. Calculations from Attaccalite et al 
(2002) are shown as crossed circles. Adopted from Zhu et al (2003). 



one can obtain the magnetization M by integrating the induced current over n s . The 
magnetic susceptibility can then be determined from the slope of the M(B) versus 
B dependence at small fields. To date, Prus et al have reported data for one sample 
which becomes insulating at a relatively high electron density (1.3 • 10 11 cm' 2 ), and 
the data obtained below this density are irrelevant for the metallic regime we are 
interested in. For this reason, the data collected so far do not provide information 
about the most interesting regime just above n x . 

In GaAs/AlGaAs heterostructures, a similar strong increase of the spin 
susceptibility at ultra-low carrier densities has now been established based on an 
analysis of the Shubnikov-de Haas oscillations (Zhu et al 2003). The results are 
shown in Fig. [23 by solid symbols; x increases by more than a factor of two as the 
density decreases. Zhu et al suggested an empirical equation \ n 7°' 4 to describe 
their experimental data, which works well in the entire range of electron densities 
spanned. Although \ tends toward a divergence, it is not clear from these experiments 
whether it does so at a finite density. We note that due to the lower effective mass, 
higher dielectric constant and the absence of valley degeneracy, the ratio r* between 
Coulomb and Fermi energies in GaAs/AlGaAs is an order of magnitude smaller than 




Figure 23. Left hand panel: magnetoresistance of GaAs/AlGaAs as a function of 
in parallel-field for different electron densities indicated in unit of 10 10 cm -2 ; the 
positions of the magnetic fields, B c , required for full polarization are indicated 
following Tutuc et al (2002). Right hand panel: B c as a function of density. 
Adopted from Zhu et al (2003). 



in silicon MOSFETs at the same electron density; therefore, to reach the same relative 
interaction strength, samples are required which remain metallic at densities about 
two orders of magnitude lower than in silicon, i.e., < 10 9 cm~ 2 . The currently 
accessible electron densities in GaAs/AlGaAs heterostructures are still too high to 
reliably establish the limiting behaviour of x- 

The open circles in Fig. 1221 denote x( n s) derived from an alternative method for 
measuring B c based on a determination of the parallel magnetic field corresponding 
to the "knee" of the magnetoresistance curves, as shown in the left hand panel of 
Fig. 1231 Earlier, this method yielded an anomalous and quite puzzling decrease of the 
susceptibility with decreasing density in both hole (Tutuc et al 2002) and electron 
(Noh et al 2002) systems in GaAs/AlGaAs. This was in disagreement with findings 
in Si MOSFETs, and argued against any spontaneous spin polarization. These results 
are now believed to be in error for a number of possible reasons. The field for full 
spin polarization, marked by vertical bars in Fig. 1231 (left hand panel) and plotted as a 
function of n s in Fig. 1231 fright hand panel), decreases with decreasing carrier density 
and exhibits an apparent saturation below n s w 0.23 TO 10 cm~ 2 . Calculation based on 
this saturation field would yield a spin susceptibility that decreases below this density. 
However, the saturation field may be constant below this density due to the fact that 
the experiments are performed at a temperature which is too high; if the temperature 
were decreased further, the saturation field would probably continue to decrease, 
consistent with the decrease of B a shown in Fig. [5UI (a). Another possible source 
of error is the finite thickness of the electron layer in GaAs/AlGaAs heterostructures, 
which leads to an increase in the effective mass with increasing parallel magnetic field 
(Batke and Tu 1986). Indeed, it has recently been shown by Zhu et al (2003) that 
the conclusion that the spin susceptibility decreases with decreasing carrier density is 
erroneous and originates from the fact that the effective mass depends on magnetic 
field. The effect of the finite thickness on the spin susceptibility was studied in detail 
by Tutuc et al (2003). 




Figure 24. The parameters g*m*, m* , and g* for different silicon samples as 
a function of r s (dots). The solid line in (a) shows the data of Okamoto et al 
(1999). The solid and open dots (b) and (c) correspond to two different methods 
of finding m* (see the text). The solid and dashed lines in (b) are polynomial fits 
for the two functions for m*(r s ). Adopted from Pudalov et al (2002). 



4-2. Effective mass and g-factor 

In principle, the increase of the spin susceptibility could be due to an enhancement 
of either g* or m* (or both). To obtain g* and m* separately, Pudalov et al (2002b) 
performed measurements of SdH oscillations in superimposed and independently 
controlled parallel and perpendicular magnetic fields. Their results are shown in 
Fig. 12] The data were taken at relatively high temperatures (above 300 mK), where 
the low-n s resistance depends strongly on temperature. Therefore, the conventional 
procedure of calculating the effective mass from the temperature dependence of 
the amplitude of the SdH oscillations is unreliable, as it assumes a temperature- 
independent Dingle temperature. Pudalov et al considered two limiting cases: a 
temperature-independent To, and a Tq that linearly increases with temperature; two 
sets of data based on these assumptions are plotted in Fig. (b) and (c). Within 
the uncertainty associated with the use of these two methods, both g* and m* were 
found to increase with decreasing n s . 

Shashkin et al (2002) used an alternative method to obtain g* and m* separately. 
They analyzed the data for the temperature dependence of the conductivity in zero 
magnetic field using the recent theory of Zala et al (2001). According to this theory, 
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Figure 25. For a low-disordered silicon MOSFET, the temperature dependence 
of the normalized conductivity at different electron densities (indicated in units of 
10 11 cm -2 ) above the critical electron density for the metal-insulator transition. 
The dashed lines are fits to the linear portions of the curves. The inset shows 
1/A*(n s ) (see Eq.[3J and B c (n a ) (open and solid circles, respectively). The dashed 
lines are continuations of the linear fits, which extrapolate to the critical electron 
density for the metal-insulator transition. From Shashkin et al (2002). 



a is a linear function of temperature: 

= 1 - A*k B T, (3) 

where the slope, A*, is determined by the interaction-related parameters: the Fermi 
liquid constants Ftf and Ff: 

(1 + 8F$)g*m* g* _ 1 

irTi 2 n s ' go 1 + F§ ' m b 



A- = - ^::° » '" , y ~ = T ^r a , ^ = l + Fl (4) 



(here go — 2 is the "bare" g-factor.) These relations allow a determination of the 
many-body enhanced g* factor and mass m* separately using the data for the slope 
A* and the product g*m*. 

The dependence of the normalized conductivity on temperature, a(T)/ao, is 
displayed in Fig. at different electron densities above n c ; the value of cto used 
to normalize a was obtained by extrapolating the linear interval of the o~(T) curve to 
T = 0. The inverse slope 1/A* (open circles) and B c (n s ) (solid circles) are plotted as 
a function of n s in the inset to Fig. [53 Over a wide range of electron densities, 1/A* 
and /ibB c are close to each other; the low-density data for 1/A* are approximated 
well by a linear dependence which extrapolates to the critical electron density n c in a 
way similar to B c . 

Values of g* /go and m* /mb determined from this analysis are shown as a function 
of the electron density in Fig. |2U In the high n s region (relatively weak interactions) , 
the enhancement of both g and m is relatively small, both values increasing slightly 
with decreasing electron density in agreement with earlier data (Ando et al 1982). 
Also, the renormalization of the g-factor is dominant compared to that of the effective 
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mass, consistent with theoretical studies (Iwamoto 1991; Kwon et al 1994; Chen and 
Raikh 1999). In contrast, the renormalization at low n s (near the critical region), 
where r s 3> 1, is much more striking. As the electron density is decreased, the 
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Figure 26. The effective mass (circles) and g factor (squares) in a silicon 
MOSFET, determined from the analysis of the parallel field magnetoresistance 
and temperature-dependent conductivity, versus electron density. The dashed 
lines are guides to the eye. From Shashkin et al (2002). 
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Figure 27. For a silicon MOSFET, the effective mass obtained from an analysis 
of p(B||) (dotted line) and from an analysis of SdH oscillations (circles). Data 
of Shashkin et al (2003a). The upper x-axis shows r* calculated using the latter 
value for the effective mass. The shaded areas correspond to metallic regimes 
studied in papers by Shashkin et al (2001b) and Prus et al (2002), as labelled. 
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Figure 28. For silicon MOSFETs, the effective mass versus the degree of spin 
polarization for the following electron densities in units of 10 11 cm -2 : 1.32 (dots), 
1.47 (squares), 2.07 (diamonds), and 2.67 (triangles). The dashed lines are guides 
to the eye. From Shashkin et al (2003a). 



renormalization of the effective mass increases markedly with decreasing density while 
the g factor remains relatively constant. Hence, this analysis indicates that it is the 
effective mass, rather than the g-factor, that is responsible for the strongly enhanced 
spin susceptibility near the metal-insulator transition. To verify this conclusion, 
Shashkin et al (2003a, 2003b) used an independent method to determine the effective 
mass through an analysis of the temperature dependence of the SdH oscillations 
similar to the analysis done by Smith and Stiles (1972) and Pudalov et al (2002b), 
but extended to much lower temperatures where the Dingle temperature is expected 
to be constant and the analysis reliable. In Fig. |23 the effective mass obtained by 
this method is compared with the results obtained by the procedure described above 
(the dotted line). The quantitative agreement between the results obtained by two 
independent methods supports the validity of both. The data for the effective mass are 
also consistent with data for spin and cyclotron gaps obtained by magnetocapacitance 
spectroscopy (Khrapai et al 2003). 

To probe a possible connection between the effective mass enhancement and spin 
and exchange effects, Shashkin et al (2003a, 2003b) introduced a parallel magnetic 
field component to align the electrons' spins. In Fig.[2Hl the effective mass is shown as 
a function of the spin polarization, P = {B\_ + B?, ) 1 / 2 /i? c . Within the experimental 
accuracy, the effective mass does not depend on P. Therefore, the enhancement of 
m* near the MIT is robust, and the origin of the mass enhancement has no relation 
to the electrons' spins and exchange effects. 

4-3. Electron- electron interactions near the transition 

The interaction parameter, r* = Ec/Ep, is generally assumed to be equal to the 
dimensionless Wigner-Seitz radius, r s = l/(7rn s ) 1 / 2 as, and, hence, proportional to 
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Figure 29. p(T) in a moderately disordered silicon MOSFET. Densities from 
top to bottom are n s = 1.20 to 1.44 (in 0.03 steps), 1.56, 1.8 and 2.04- 10 11 cm" 2 . 
From Prus et al (2002). 

— 1/2 

n s . However, this is true only if the effective mass does not depend on n s ; close to 
the transition, where the effective mass is strongly enhanced, the deviations from the 
square root law become important, and the interaction parameter is larger than the 
Wigner-Seitz radius by a factor of m* /m& (or 2m* jm^ in silicon MOSFETs, where an 
additional factor of 2 derives from the valley degeneracy). 

Near the metal-insulator transition, r* grows rapidly due to the sharp increase of 
the effective mass, as shown on the upper ir-axis in fjgurc l2"YI here r* is calculated using 
the effective mass obtained from the analysis of the SdH oscillations. We note that 
since at low n s this method yields a somewhat smaller effective mass than the method 
based on the analysis of p(B\\), the plotted values represent the most conservative 
estimate for r* . Due to the rapid increase of m* near the transition, even small changes 
in the electron density may lead to large changes in r* . For example, the transition to 
a localized state in a low-disordered sample used by Shashkin et al (2001b) occurs at 
n s = 0.795- 10 11 cm -2 (see Fig. 1151) . while in a sample of lesser quality used by Prus et 
al (2002), it occurs at 1.44 • 10 11 cm~ 2 (see Fig. I29p . The areas corresponding to 
the metallic regimes studied in these papers are shaded in Fig. [23 the corresponding- 
values of r* differ dramatically (the actual values of r* in the study by Shashkin et al 
are not known because the effective mass has not been determined near the transition; 
the lower boundary for r* is 50). The difference in r* may account for the qualitative 
change in the behaviour of the resistance in the two samples: in the more disordered 
sample, there is no temperature- independent curve (the separatrix), and some of the 
curves which look "metallic" at higher temperature are insulating below a few hundred 
mK, suggesting that localization becomes dominant in this sample as T — > 0. 

5. HOW DOES ALL THIS FIT THEORY? 

The possibility that a B = metallic state in 2D can be stabilized by interactions 
was first suggested by Finkelstein (1983, 1984) and Castellani et al (1984). In this 
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theory, the combined effects of disorder and interactions were studied by perturbative 
renormalization group methods. It was found that as the temperature is decreased, 
the resistivity increases and then decreases at lower temperatures, suggesting that the 
system is approaching a low-temperature metallic state. An external magnetic field, 
via Zeeman splitting, drives the system back to the insulating state. These predictions 
of the theory are in qualitative agreement with experiments. 

However, an interaction parameter scales to infinitely large values before zero 
temperature is reached, and the theory thus becomes uncontrolled; this scenario has 
consequently not received general acceptance. It should also be noted that the theory 
may not be applicable to the current experiments since it was developed for the 
diffusive regime: T <C h/r, where r is the elastic mean- free time extracted from the 
Drude conductivity (Boltzmann constant is assumed to be equal to 1 throughout this 
section). This condition corresponds to the low-temperature limit T <C Tpp/(h/e 2 ). 
Since the Fermi temperature, Tp, is rather low at the small carrier densities considered 
here, the above condition is satisfied only close to the transition, where p becomes 
of the order of h/e 2 . In the experiments, however, the characteristic decrease of 
the resistance with decreasing temperature often persists into the relatively high- 
temperature ballistic regime T > H/t (or T > T F p/(h/e 2 )). Altshuler et al (2001) 
and Brunthalcr et al (2001) have interpreted this observation as evidence that the 
mechanism responsible for the strong temperature dependence cannot originate from 
quantum interference. 

For the ballistic region, calculations in the random phase approximation were 
carried out two decades ago by Stern (1980), Gold and Dolgopolov (1986), Das Sarma 
(1986) and were recently refined by Das Sarma and Hwang (1999). These theories 
predict a linear temperature dependence for the conductivity with metallic-like slope 
(da/dT < 0) regardless of the strength of the interactions. The spin degree of 
freedom is not important in this theory and enters only through the Fermi energy, 
which is a factor of two larger for the spin-polarized than for the unpolarized system. 
Therefore, the application of a (parallel) magnetic field does not eliminate the metallic 
temperature dependence. 

The two limits — diffusive and ballistic — had until recently been assumed to 
be governed by different physical processes. However, Zala et al (2001) have shown 
that the temperature dependence of the conductivity in the ballistic regime originates 
from the same physical process as the Altshuler-Aronov-Lee corrections: coherent 
scattering of electrons by Friedel oscillations. In this regime the correction is linear 
in temperature, as is the case for the results mentioned in the previous paragraph. 
However, the value and even the sign of the slope depends on the strength of electron- 
electron interaction, the slope being directly related to the renormalization of the 
spin susceptibility (see also Gold 2001, 2003). By aligning the spins, a magnetic field 
causes a positive magnetoresistance and changes the temperature dependence of the 
conductivity from metallic-like to insulating-like (Herbut 2001; Zala et al 2002; Gold 
2003), in agreement with experiments. 

5.1. The diffusive regime: Renormalization group analysis 

It is well known that in the diffusive limit (TV -C h), electron-electron interactions in 
two dimensional disordered systems lead to logarithmically divergent corrections to 
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the conductivity given by: 



6a = -^n ln (^ 



1 + 3 1 



In(l + F ») 



(•5) 



where Ffi is the interaction constant in the triplet channel which depends on the 
interaction strength. The sign of this logarithmically divergent correction may be 
positive (metallic-like) or negative (insulating- like), depending on the value of Fq. 

Experimentally, the diffusive regime is realized in a relatively narrow range of 
electron densities near the metal-insulator transition, where the resistivity is of the 
order of h/e 2 . Punnoose and Finkelstein (2002) have convincingly demonstrated 
that in this region, the temperature dependence of the resistivity can be understood 
within the renormalization group theory describing the effect of the electron-electron 
interactions on the propagation of diffusive collective modes, with the delocalizing 
component becoming dominant in dilute 2D systems. 

In 2D, the renormalization group equation describing the evolution of the 
resistivity has been derived previously by Finkelstein: 



dg 2 

~J7 = 9 



1 + 1-3 ( i±2*ln(l+ 72 )-l 

72 



(6) 



Here £ = — IniTr/h), the dimensionless parameter g is the conductance in units 
of e 2 /irh and 72 is the coupling constant. The first term in the square brackets 
of Eq. El corresponds to the weak localization correction (quantum interference; 
Abrahams et al 1979), while the second term is the contribution of the singlet density 
mode which is due to the long range nature of the Coulomb interaction (Altshuler, 
Aronov and Lee 1980). The last term describes the contribution of the three triplet 
modes. Note that the last two terms have opposite signs, favouring localization and 
derealization, respectively. The resulting dependence <?(£) becomes delocalizing when 
72 > 7 2 = 2.04. This requires the presence of rather strong electron correlations. 

In the case of two distinct valleys (as in (100) silicon MOSFETs), Eq. Elcan be 
easily generalized to: 



dg 2 
-77 = .9 

d£ 



1 _ 15 ( l±^ln(l + 72 )-l 

72 



(7) 



The difference between the numerical factors in Eq. and Eq. [7] results from the 
number of degrees of freedom in each case. The weak localization term becomes twice 
as large. The difference in the number of the multiplet modes increases the coefficient 
of the 72 term from 3 to 15. As a result of these modifications, the value of 72 required 
for the dependence <?(£) to become delocalizing is reduced to 72 = 0.45. which makes 
it easier in the case of two valleys to reach the condition where the resistivity decreases 
with decreasing temperature. 

In conventional conductors the initial values of 72 are small, and the net effect 
is in favour of localization. In dilute systems, however, this value is enhanced due to 
electron-electron correlations. In addition, 72 also experiences logarithmic corrections 
due to the disorder (Finkelstein 1983, 1984; Castellani et al 1984, 1998). The equation 
describing the renormalization group evolution of 72 is the same for the case of one 
and two valleys: 

d l2 _ (I+72) 2 , * 

It follows from this equation that 72 increases monotonically as the temperature is 
decreased. Once it exceeds 7|, the resistivity passes through a maximum. Even though 



Metal-insulator transition in 2D 



38 




Figure 30. Data for silicon MOSFETs at n s = 0.83, 0.88 and 0.94 ■ 10 11 cm -2 
from Pudalov et al (1998) are scaled according to Eq.|2] The solid line corresponds 
to the solution of the renormalization group, Eqs. |7| and [S] no adjustable 
parameters were used in this fit. From Punnoose and Finkelstein (2002). 



the initial values of g and 72 are not universal, the flow of g can be described by a 
universal function R(rj) (Finkelstein, 1983): 

9 = 9ma X -R(l) and 77 = g max ln(T max /T) , (9) 

where T max is the temperature at which g reaches its maximum value <? max . 

In Fig. |20l the theoretical calculations are compared with the experimental data 
obtained by Pudalov et al (1998). Since the renormalization group equations were 
derived in the lowest order in g and cannot be applied in the critical region where 
g > 1, only curves with maximum g ranging from g max ~ 0.3 to g ma , x ~ 0.6 are shown 
in the figure. The decrease in the resistivity by up to a factor of five, as well as 
its saturation, are both captured in the correct temperature interval by this analysis 
without any adjustable parameters. 

In agreement with experiment, this universal behaviour can be observed only in 
ultra clean samples (with negligible inter-valley scattering) and will not be found in 
samples that are only moderately clean. In disordered silicon MOSFETs, a description 
in terms of an effective single valley is relevant, and the large value of 7J = 2.04 makes 
it difficult for the non-monotonic p{T) to be observed as the initial values of 72 are 
usually much less than 2.04. As a result, for g near the critical region, the resistivity 
becomes insulating-like instead of going through the maximum. 

Therefore, in ultra clean silicon MOSFETs, the behaviour of the resistivity not 
far from the transition is quantitatively well described by the renormalization group 
analysis that considers the interplay of the electron-electron interactions and disorder 
when the electron band has two distinct valleys. The theory in this case remains in 
control down to exponentially low temperatures, unlike the single valley case, where 
72 diverges at r\ w 1, or at temperatures just below T max . A parallel magnetic 
field provides a Zeeman splitting and gives rise to positive magnetoresistance (Lee 
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Figure 31. Total interaction correction to the conductivity. The curve for 
Fq = corresponds to the universal behaviour of a completely spin polarized 
electron gas. From Zala et al (2001). 



and Ramakrishnan 1982; Finkclstcin 1984; Castcllani et al 1998). In a very strong 
magnetic field, when the electrons are completely spin-polarized, the system becomes 
"spinless". In this case, the system always scales to an insulator, and in the weak 
disorder limit, a universal logarithmic temperature dependence is expected: 

a{T) = cr + {e 2 /Trh){2 - 2 In 2) \u{Tt/K). 



5. 2. Farther from the transition ( the ballistic regime ) 

Far from the transition {n s n c ), the inverse scattering time r _1 is often smaller than 
the temperature at which the experiments are performed, and one is in the ballistic 
regime: Tr 3> h. The interaction theory by Zala et al (2001) considers coherent 
electron scattering off the Friedel oscillations. The phase of the Friedel oscillation 
is such that the wave scattered from the impurity interferes constructively with the 
wave scattered from the oscillation, leading to a quantum correction to the Drude 
conductivity, uq. As has already been mentioned, this correction (Act), which is 
logarithmic in the diffusive regime, becomes linear in the ballistic regime: 

a. ( t) = 4^ ( i + j*u = cto ( i + ^ do 

where Fq is the Fermi liquid interaction parameter in the triplet channel. As in the 
diffusing regime, the sign of dp/dT depends on the constant Fq . 

The total correction to the conductivity is shown in Fig. 1311 for different values of 
Fq . The divergence at low temperature is due to the usual logarithmic correction 
(Altshuler, Aronov and Lee 1980; Finkelstein 1983, 1984; Castellani et al 1984). 
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Figure 32. (a) Ap as a function of the dimensionless temperature (T/Tp) at 
different p in a p-type GaAs/AlGaAs heterostructure. (For clarity, curves in (a) 
and (b) arc offset vertically.) (b) The same data as in (a) plotted as conductivity, 
with linear fits, (c) Fermi liquid parameter versus hole density. Open symbols 
show the result obtained from the analysis of p(T) at zero magnetic field; closed 
symbols show results obtained from the parallel-field magnetoresistance. From 
Proskuryakov et al (2002). 



Although the exact value of Ffi cannot be calculated theoretically (in particular, 
its relation to the conventional measure of the interaction strength, r s , is unknown 
for r s > 1), it can in principle be found from a measurement of the Pauli spin 
susceptibility. The correction to the conductivity is almost always monotonic, except 
for a narrow region —0.45 < Ffi < —0.25. 

As in the diffusive limit, magnetic field in the ballistic regime suppresses the 
correction in the triplet channel in Eq. (|10|l . resulting in a universal, positive correction 
to the Drude conductivity in magnetic field, <Jq , and hence in the insulating-like 
behaviour of cr(T): 

Aa = o-Z^L &t B > B s , (11) 

where B s is the field corresponding to full spin polarization of the 2D system. 

Proskuryakov et al (2002) were the first to perform a quantitative comparison 
of experimental data (obtained on p-type GaAs/AlGaAs) with the above theory. 
Their data for the temperature-induced corrections to the conductivity are plotted 
in Fig. |22(b). (To extract corrections to the conductivity, Proskuryakov et al used 
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Figure 33. For parameters A = 1 K and Fq = —0.15, the magneto-conductivity 
of a silicon MOSFET: (a) For different temperatures at a fixed density n s = 
2.47 ■ 10 11 cm" 2 ; T = 0.25, 0.5, 0.8, 1.0, 1.3, 1.6, 2.6 and 3.6 K (from top), 
(b) For different densities at a fixed temperature of 1 K; n s = 3.3, 3.0, 2.75, 2.47, 
2.19, 1.92 and 1.64 • 10 11 cm -2 (from top), (c) The conductivity as a function of 
temperature for the densities n s = 3.3, 3.0, 2.47, 2.19, 1.92 and 1.64 ■ 10 11 cm -2 
(from top). Symbols and solid lines denote data and theory, respectively. From 
Vitkalov et al (2003). 

the Ap(T) dependence shown in Fig. [221 (a) obtained from raw resistivity data by 
subtracting the contribution of phonon scattering.) A linear fit of A<t(T) yields the 
parameter Fq shown in Fig. [221 (c) for different p. The interaction constant obtained 
is negative and provides the metallic slope of <j{T). As expected, the slope increases 
with decreasing density. When extrapolated to the crossover point from metallic to 
insulating behaviour (p « 1.5 • 10 10 cm" 2 ), the interaction constant is much higher 
than the value of Fq = — 1 for the Stoner instability. The transition from metallic- like 
to insulating-like o~(T) in their sample is unlikely to be related to any dramatic change 
in magnetic properties and is probably caused by Anderson localization. 
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Figure 34. Schematic phase diagram corresponding to the numerical results 
of Attaccalite et al 2002. The symbols WC, FFL and PFL correspond to the 
Wigner crystal, ferromagnetic Fermi liquid and paramagnetic Fermi liquid phases, 
respectively. 



Comparison between the experimentally measured <r(T) in the ballistic regime 
with the predictions of Zala et al was also carried out in silicon MOSFETs by 
Shashkin et al (2002; see sec.lOJl. Vitkalov et al (2003) and Pudalov et al (2003) 
(see also Das Sarma and Hwang 2003 and Shashkin et al 2003c) and in p-type SiGc 
heterostructures by Coleridge et al (2002). As shown in Fig. 123(a) and (b), Vitkalov et 
al reported quantitative agreement between the theory and their magnetoresistance 
data in silicon MOSFETs. However, the values of the valley splitting A and Fermi 
liquid parameter Fq required to obtain fits to the field dependence yield curves that 
are inconsistent with the observed temperature dependence in zero field, as shown in 
Fig. [221 (c). This was attributed to the neglect of additional scattering terms that 
affect the temperature dependence more strongly than the field dependence. Despite 
this quantitative discrepancy, the theory of Zala et al provides a reasonable description 
of the conductivity of strongly interacting 2D systems in the ballistic regime. 

5.3. Approaches not based on Fermi liquid 

The calculations of Finkelstein (1983, 1984), Castellani et al (1984) and Zala et al 
(2001) use the Fermi liquid as a starting point. As discussed earlier, r s becomes so 
large that the theory's applicability is in question near the transition (Chamon et al 
2001; for a review, see Varma et al 2002). Moreover, the behaviour of the effective 
mass reported by Shashkin et al (2003a, 2003b) in the vicinity of the transition is 
unlikely to be consistent with a Fermi liquid model. 

Very little theory has been developed for strongly interacting systems for which 
r s is large but below the expected Wigner crystallization. Several candidates have 
been suggested for the ground state of the strongly interacting 2D system, among 
them (i) a Wigner crystal characterized by spatial and spin ordering (Wigner 1934), 

(ii) an itinerant ferromagnet with spontaneous spin ordering (Stoner 1946), and 

(iii) a paramagnetic Fermi liquid (Landau 1957). According to numerical simulations 
(Tanatar and Ceperley 1989), Wigner crystallization is expected in a very dilute 
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Figure 35. Phase diagram of the 2D electron system at T = suggested by 
Spivak (2002). The symbols WC and FL correspond to the Wigner crystal and 
the Fermi liquid phases, respectively. The shaded regions correspond to phases 
which are more complicated than the bubble and the stripe phases. 



regime, when r s reaches approximately 35. This value has already been exceeded 
in the best p-type GaAs/AlGaAs heterostructures; however, no dramatic change 
in transport properties has occurred at the corresponding density. Recent detailed 
numerical simulations (Attaccalite et al 2002) have predicted that in the range of 
the interaction parameter 25 < r s < 35 prior to the crystallization, the ground 
state of the system becomes an itinerant ferromagnet. The corresponding schematic 
phase diagram is shown in Fig. |MJ As discussed earlier, there are experimental 
indications that a spontaneous spin polarization may occur at a finite electron density 
in silicon MOSFETs (and possibly in GaAs/AlGaAs heterostructures). Chakravarty et 
al (1999) have proposed a more complicated phase diagram where the low-density 
insulating state is a Wigner glass, a phase that has quasi-long-range translational 
order and competing ferromagnetic and antiferromagnetic spin-exchange interactions. 
The transition between insulating and metallic states within this theory is the melting 
of the Wigner glass, where the transition is second order due to the disorder. 

Spivak (2002) predicted the existence of an intermediate phase between the Fermi 
liquid and the Wigner crystal with a first order transition in a clean electron system. 
The suggested phase diagram is shown in Fig. ESI I n analogy with He 3 , where m* 
increases approaching the crystallization point, Spivak (2001) suggested in an earlier 
paper that the renormalization of m* is dominant compared to that of the g-factor 
as the transition is approached, and that m* should increase with magnetic field. 
Although the increase of the mass is in agreement with the experimental results of 
Shashkin et al (2002, 2003a, 2003b), the suggested increase of m* with the degree 
of spin polarization is not confirmed by their data. The strong increase of the 
effective mass near crystallization also follows from the approach used by Dolgopolov 
(2002), who has applied Gutzwiller's variational method (Brinkman and Rice 1970) to 
silicon MOSFETs, and from the dynamical mean-field theory (Tanaskovic et al 2003). 
However, the predicted dependence of m* on the degree of spin polarization is again 
at odds with the experiment. Dharma-wardana (2003) has argued that in two-valley 
systems, correlation effects outweigh exchange, and a coupled- valley state is formed 
at r s k 5.4 leading to strong enhancement of the effective mass, which in this case is 
only weakly dependent on the degree of spin polarization. 

6. CONCLUSIONS 

Significant progress has been achieved during the past few years in understanding 
the metallic state found a decade ago in strongly interacting 2D systems. There 
is now considerable evidence that the strong metallic temperature dependence of 
the resistance in these systems is due to the delocalizing effects of strong electron- 
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electron interactions. In ultra-clean silicon systems, the temperature dependence of 
the resistivity is universal near the metal-insulator transition and is quantitatively 
well described by the renormalization group theory; deep in the metallic state, in 
the ballistic regime, the temperature dependence of the resistance can be explained 
by coherent scattering of electrons by Friedel oscillations. In both cases, an external 
magnetic field quenches the delocalizing effect of interactions by aligning the spins, 
causing a giant positive magnetoresistance. 

The metal-insulator transition is not yet understood theoretically. In silicon 
MOSFETs, various experimental methods provide evidence for a sharp increase 
and possible divergence of the spin susceptibility at some finite sample-independent 
electron density, n x , which is at or very near the critical density for the MIT in 
high mobility samples. Unlike the Stoner instability which entails an increase in the 
g-factor, the increase in the susceptibility in these systems is due to an increase of 
the effective mass. The effective mass is, in turn, found to be independent of the 
degree of spin polarization, implying that the increase is not due to spin exchange, in 
disagreement with the Fermi liquid model. A similar increase of the spin susceptibility 
is observed in GaAs/AlGaAs heterostructures, but it is not yet clear whether or not 
it points to a spontaneous spin polarization at a finite carrier density. 

The fact that the B = metal- insulator transition in the least disordered silicon 
samples occurs at or very close to n x indicates that the transition in such samples 
is a property of a clean electron system and is not driven by disorder. Quantum 
localization appears to be suppressed near the transition in these systems. In lower 
mobility samples, the localization transition occurs at electron densities much higher 
than n x and may be driven by disorder. 

In closing, we note that most of the work done in these dilute two-dimensional 
systems concerns the transport behaviour, as transport measurements are relatively 
straightforward (despite problems associated with high impedance electrical contacts 
at low densities, the need to ensure proper cooling of the electrons system and so 
on). Studies done to date, many of which are reviewed here, include measurements 
of the resistivity as a function of temperature and magnetic field, the Hall coefficient, 
Shubnikov-de Haas oscillations and measurements of noise. Results have also been 
reported for the compressibility (Dultz and Jiang 2000; Hani et al 2000, 2001; see 
also Si and Varma 1998 and Fogler 2003) and capacitive measurements (Khrapai et 
al 2003) from which one can obtain information about the chemical potential, the 
density of states, and which provide a measure of the magnetization, as discussed in 
sec. 14.11 

Many properties that would yield crucial information have not been investigated. 
Thermodynamic measurements such as specific heat and direct measurements of 
magnetization would be particularly illuminating; however, these are very difficult, 
if not impossible, to perform at this time due to the very small number of 
electrons available in a dilute, two-dimensional layer. Other experiments which could 
provide valuable information include tunnelling and different resonance techniques. 
Measurements of one, several, or perhaps all of these may be required for a full 
understanding of the enigmatic and very interesting behaviour of strongly interacting 
electrons (or holes) in two dimensions. 
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